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Full Jet Reconstruction in Heavy Ion Collisions: Prospects and Perils 
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Full jet reconstruction in heavy ion events has been thought to be difficult due to large multiplicity backgrounds. 
A new generation of jet reconstruction algorithms to search for new physics in high luminosity p+p collisions 
at the LHC is developed to precisely measure jets over large backgrounds caused by pile up. From simulations 
it turns out, this new generation of reconstruction algorithms are also applicable in the heavy ion environment. 
We review the latest results on jet-medium interactions as seen in A+A collisions at RHIC, focusing on the new 
techniques for full jet reconstruction. 



1. Introduction 

Jets must be well-defined, measurable from the 
hadronic final-states and calculable in perturbative 
QCD (pQCD) from the partonic states 0, Q. The 
precise measurements of inclusive jet cross sections 
are performed at many hadronic and leptonic colliders 
to check in detail the pQCD calculations, to help de- 
termine parton distribution functions and to look for 
new physics. These measurements are found to be in a 
very good agreement with the different pQCD calcula- 
tions using various parton distribution functions. The 
robustness of the theoretical calculations on jet cross 
sections in p + p collisions motivates the use of jets as 
direct probes of partonic energy loss in the hot QCD 
matter generated in ultra-relativistic heavy ion colli- 
sions at RHIC and in near future at the LHC through 
their interaction and energy loss in the medium ("jet 
quenching") 

Until recently, inclusive hadron distributions and 
di-hadron correlations at high transverse momentum 
were utilized to measure jet quenching at RHIC indi- 
rectly to avoid the complex backgrounds of high mul- 
tiplicity heavy ion events. However, these measure- 
ments selects jet fragmentation particles that are bi- 
ased towards the population of jets that has the least 
interaction with the medium. Various models with a 
wide range of parameters are able to describe these 
measurements however they are leading to only lim- 
ited constraints upon the underlying physics [J, [5|, [6|]. 
Full jet reconstruction in A+A collisions can overcome 
theses geometric biases as the energy flow is measured 
independently of the fragmentation details. Jet recon- 
struction at the partonic level with significantly re- 
duced biases enables this study, with qualitatively new 
observables such as energy flow, jet substructure and 
fragmentation functions that can be measured in mul- 
tiple channels (inclusive, di-jets, h-jets and gamma- 
jets). 

The detector upgrades together with the increased 



beam luminosities of RHIC and data recording capa- 
bilities of PHENIX and STAR, enable jet reconstruc- 
tion in heavy ion collisions for the first time [7|. In 
this article we review the prospects and the perils of 
the new techniques developed for full jet reconstruc- 
tion in heavy ion collisions at RHIC and the latest 
results obtained with these new tools. The experimen- 
tal details of jet reconstruction measurements utilizing 
the STAR and PHENIX experiments can be found in 
@, i, ES EH for the inclusive spectra, [H, El for the 
underlying event, and 0, EH for the accompanying 
jet fragmentation studies in heavy ion collisions. 



2. Prospects 

Various jet reconstruction algorithms have been de- 
veloped for both leptonic and hadronic colliders dur- 
ing the last 20 years. Here we will briefly discuss the 
cone and sequential recombination algorithms that are 
used for the STAR analysis. Further algorithmic de- 
tails of cone, sequential recombination and Gaussian 
filtering can be found in [H El El, El, HI EH, S3 
and the references therein. 

2.1. Jet Reconstruction Algorithms 

The cone algorithm is based on the simple picture 
that a jet consists of a large amount of hadronic energy 
in a small angular region. Therefore, the main method 
for the cone algorithm is to combine particles in 7/ — <j) 
space with their neighbors within a cone of radius R 
(R = y/ A(j) 2 + An 2 ). To accommodate higher-order 
processes and to optimize the search and effectiveness 
of jet finding, splitting, merging, and iteration steps 
can be used. 

Unlike the cone algorithm, the sequential recombi- 
nation algorithms combine pairs of objects relative to 
the closeness of their pt- Particles are merged into 
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a new cluster via successive pair-wise recombination. 
While the lowest pr particles are the starting point 
for clustering in the fc-r algorithm, in the anti-fcr al- 
gorithm, recombination starts with the highest mo- 
mentum particles. Due to these different approaches, 
in the kr algorithm, arbitrarily shaped jets are al- 
lowed to follow the energy flow, resulting in less bias 
on the reconstructed jet shape than with the anti-fc^ 
or cone algorithm which are more or less restricted to 
a circular shape [Til ]. 

Motivated by the need for precision jet measure- 
ments in the search for new physics in high lumi- 
nosity p+p collisions at the LHC, a new approach 
to jet reconstruction and background subtraction was 
developed [23j . A key feature in this procedure 
is a new QCD inspired algorithm for separating jets 
from the large backgrounds due to pile up. As it 
turns out from simulations, these improved techniques 
can also be used in heavy ion environments where 
the background subtraction is essential for jet mea- 
surements. Sequential recombination algorithms (kr, 
anti-kx and Cambridge / Aachen ( C AMB) ) encoded in 
the Fast, Jet suite of programs 0, HI]) along with a 
seedless infrared-safe cone algorithm (SISCone) [25j ] 
are utilized to search for jets in p+p [26| and Au+Au 
collisions collected by the STAR experiment. 

An alternative seeded cone algorithm was also ex- 
plored previously by STAR experiment for Au+Au 
collisions in order to avoid instabilities in cone-finding 
due to large heavy ion background [H, [3] • In addition 
to recombination and cone algorithms, the Gaussian 
filtering type of algorithm that simply extracts jets as 
local maxima in the f] — (f) space by linearly filtering 
particles is also used by the PHEN1X experiment to 
extract jets successfully in p+p and Cu+Cu collisions 



3. Perlis 



The complex heavy ion background makes the full 
jet reconstruction a challenging task. Here we dis- 
cuss the fundamental assumptions required by the full 
jet reconstruction in these complex environments and 
the new biases that can be introduced. These per- 
ils require further investigation of the algorithmic re- 
sponses of full jet reconstruction in heavy ion collisions 
by utilizing various reconstruction algorithms and the 
implications on jet quenching in heavy ion collisions. 



3.1. Underlying Heavy Ion Background 



A fundamental assumption that the signal and the 
uniform background are two separable components is 
required by the full jet reconstruction in heavy ion col- 
lisions. However, this assumption can be violated by 
the presence of jets and the effect of these jets on the 
background estimation. For example the initial state 
radiation, even though expected to be small compared 
to jet energy, might be different in A+A relative to 
p+p. Initial state processes resulting in the enhance- 
ment of the multiplicity of the underlying events ap- 
pears to be distributed uniformly in the simpler p+p 
case [IH and therefore it is fully accounted for the esti- 
mation of the background under jets [l6j |. However in 
the Au+Au collisions, "the p+p correspondent" un- 
derlying event may be modified, possibly generating 
non-uniform structures. These non-uniformities in the 
background might be even larger due to the final state 
processes. Energy loss of the jet in matter might mod- 
ify the event shape, resulting in non-uniform struc- 
tures such as the ridge. Azimuthal and longitudinal 
anisotropy of heavy ion events will also result into 
non-uniform backgrounds. Some of these sources of 
correlated backgrounds can be brought under quanti- 
tative control by using different collision systems. On 
the other hand, other observed effects might help us 
to understand details of the jet interactions with the 
heavy ion environment and may give further insight 
into the structures that are observed in di-hadron cor- 
relations and their origins. 

With the assumption that the signal and the back- 
ground are two separable components, the background 
correction can be estimated by following three simple 
steps. The first step is measuring the jet area for 
the infrared safe algorithms as encoded in the Fast- 
Jet Suite of algorithms. An active area of each jet is 
estimated by filling an event with many very soft par- 
ticles and then counting how many are clustered into 
a given jet. The second step is measuring the diffuse 
noise (mean px per unit area in the remainder of the 
event) and noise fluctuations. Event-by-event fluctu- 
ations in the background also distort the jet spectrum 
towards larger Et due to the steeply falling depen- 
dence of the jet production on Et- This effect can be 
corrected through an unfolding procedure (i.e., decon- 
volution). The final step is correcting the jet energy 
by deconvolution of signal from the complex heavy ion 
background, using parameters that are extracted from 
measurable quantities of area, diffuse noise and noise 
fluctuations. 
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Figure 1: The pt threshold dependent comparison of inclusive jet spectra in heavy ion and p+p collisions. The filled 
triangle symbols are from MB-Trig, and filled squares are from Neinary scaled p+p collisions. In addition to pt cuts 
shown in the figures, R is set to 0.4. Q. 



3.2. Biases 

The ultimate goal of full jet reconstruction is to 
investigate the jet quenching in heavy ion collisions 
at the partonic level, without any ambiguities being 
introduced by hadronization and geometric biases of 
the inclusive spectrum and di-hadron measurements. 
However, it is possible that new biases can be intro- 
duced when reconstructing jets. For example, all jet 
algorithms have various parameters for searching and 
defining jets, and the effects of varying these param- 
eters need to be explored in detail for a full under- 
standing of jet reconstruction. 

A bias will be introduced while trying to reduce 
the effect of the background fluctuations in heavy ion 
collisions with the threshold cuts on the track mo- 
menta and calorimeter tower energies (p™* ) . Figure [T] 
shows the comparison of the jet spectra reconstructed 
by the kx algorithm with a variation of pt threshold 
cuts in Au+Au and the Nsinary scaled p+p collisions 
Q. For these results, energy resolution of the detec- 
tors and the underlying heavy ion background fluc- 
tuations were corrected with multiplicative factors of 
the jet spectrum estimated by utilizing Monte-Carlo 
model studies based on Pythia 8.107 [27|. While the 
agreement between Au+Au and Neinary scaled p+p 
jet measurements is good for the lowest value of the 
Pt cut, it is also seen to be poorer with the larger pt 
threshold cut. This suggests that the threshold cuts 
introduce biases which are not fully corrected by the 
procedures that use fragmentation models that are de- 
veloped for e + + e~ and p+p collisions. 

To enhance the recorded rate of high px particles 
and jets, events above some threshold in the electro- 
magnetic calorimeter are collected. (This threshold is 
5.4 GeV for the STAR experiment during data taking 



in years 2006 and 2007.) This is very similar to the 
case of jets that are reconstructed with seeded infrared 
unsafe algorithms. The online calorimeter tower trig- 
gers introduce a strong bias of reconstructed jets that 
are fragmenting hard in comparison to the jets that 
are reconstructed without a seed. 

The resolution parameter or cone size, which re- 
stricts the area of the jet and thereby the amount of 
energy flow, can be a harder parameter to calculate 
hence interpret in heavy ion collisions than in p+p 
collisions. If the jets are broader in the heavy ion en- 
vironment, the same resolution parameter might not 
be sufficient to recover the same fraction of jet energy 
in comparison to p+p jets. This bias needs to be in- 
vestigated by varying the resolution parameter and by 
looking into the jet profile of these jet definitions. 



4. Results 

The Figure [5] shows the comparison of the inclu- 
sive jet spectra reconstructed by kx and anti-kx se- 
quential recombination algorithms for central Au+Au 
collisions collected by the STAR experiment Sys- 
tematic uncertainties due to the unfolding procedure 
are shown as the envelopes in red and black and the 
jet energy resolution as the yellow bar. Jet spectra 
are consistent between the two sequential reconstruc- 
tion algorithms extending to 50 GeV kinematic reach. 
Corrected jet spectra reconstructed with a Gaussian 
filtering algorithm for various Cu+Cu centralities col- 
lected by the PHENIX experiment within their re- 
stricted experimental acceptance were also presented 
in this meeting [13, [TTJ] . The same algorithms that are 
used for heavy ions are also used to reconstruct jets 
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in p+p collisions. The jet spectra reconstructed in 
p+p collisions collected by PHENIX and STAR exper- 
iments (with sigma of 0.3 for gaussian filtering and the 
resolution parameter of 0.4 for the sequential recom- 
bination algorithm) all agree well with the previously 
published RHIC results using a cone al gori thm with 
split merge steps (cone radius of 0.4) [9ITT0L [28L [29l] . 



Raa above 30 GeV can be observed. This implies that 
an unbiased jet reconstruction is not reached fully. It 
is expected that for a smaller resolution parameter, 
this suppression should reach to single particle sup- 
pression at large momentum. For the case of R=0.2, 
further suppression is observed to grow to larger val- 
ues supporting this expectation [§]. 




Figure 2: Inclusive jet yield per event vs transverse jet 
energy for the central Au+Au collisions obtained by the 
sequential recombination (kx and anti-kx) algorithms |9|. 

The nuclear modification factor (Raa) for the re- 
constructed jet spectra with the resolution parameter 
of 0.4 from kx and anti-kx can be found in Figure [3l 
The envelopes shown represent the one sigma uncer- 
tainty of the deconvolution of the heavy ion back- 
ground. The total systematic uncertainty on the jet 
energy scale is around 50%, shown as the yellow bar. 
In the case of full jet reconstruction, Neinary scaling 
as calculated by a Glauber model [30( is expected if 
the reconstruction is unbiased, i.e. if the jet energy is 
recovered fully independent of the fragmentation de- 
tails, even in the presence of strong jet quenching. 
This scaling is analogous to the cross section scal- 
ing of high pt direct photon production in heavy ion 
collisions observed by the PHENIX experiment (3lj |. 
A large fraction of jets are reconstructed when using 
both kx and anti-kx sequential recombination algo- 
rithms with a resolution parameter of 0.4. Momen- 
tum dependence of the nuclear modification factor is 
also different than the observed suppression of the n 
meson Raa- However even though there are large 
systematic uncertainties, a hint of a suppression of jet 



Figure 3: Momentum dependence of the nuclear modi- 
fication factor of jet spectra reconstructed with kx and 
anti-kx algorithms (0-10% most central Au+Au divided 
by N Binary scaled p+p collisions) |J. 

The ratio of the spectra from the recoil jets of 
the di-jet coincidence measurements in 0-20% central 
Au+Au to p+p collisions is presented in Figure ffl[l5j. 
The recoil jets are selected when the triggered jets 
havepx greater than 10 GeV and 20 GeV as presented 
as the triangles and the open circles. The recoil jet 
spectra are normalized to the number of trigger jets. 
Background and trigger jet energy uncertainties are 
shown with the solid and dashed histograms. While a 
large fraction of inclusive jets are reconstructed yield- 
ing a much smaller suppression as seen in Figure [3j 
with these trigger selections biased population of re- 
coil jets are measured resulting into a suppression that 
is much more comparable to the measurement of 7r 
meson Raa- It is possible to introduce geometric bi- 
ases that are similar to the observed di-hadron mea- 
surements when reconstructed jets are selected with a 
given criteria. 

Medium effects at parton splitting can also be 
studied with fragmentation functions. Fragmentation 
function measurements are extremely challenging due 
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Figure 4: The ratio of the spectra of the recoil recoil jets 
in 0-20% central High Tower triggered Au+Au events to 
High Tower triggered p+p collisions [TEj . 



to the uncertainty in the statistical subtraction of the 
background particles and to potential biases in the re- 
construction discussed earlier. Nevertheless, the Fig- 
ure[3]is a first attempt of the z (z — pt^ adron /p^ et ) de- 
pendence of the ratio of the fragmentation functions 
from the recoil jets of the di-jet coincidence measure- 
ments in 0-20% central Au+Au to p+p collisions [H| . 
See the inset of the right panel in Figure [3] for the 
applied selection cuts for the trigger and the recoil 
jets reconstructed with the anti-fcx algorithm. Vari- 
ous uncorrelated systematic uncertainties are also in- 
cluded as the solid and the dotted histograms. Within 
the given systematic uncertainties, the ratio of the 
fragmentation functions of Au+Au to p+p show no 
significant suppression. The lack of modification in 
fragmentation functions confirms that the di-jets that 
are reconstructed in this analysis are biased towards a 
sample of unquenched jets that are only coming from 
the surface as also observed as the suppression in Fig- 
ure @] It is also possible that due to the broadening 
of the jets and the requirement of the limited resolu- 
tion parameter only part of the jet energy is recovered. 
This results in an insensitivity to expected softening 
when measuring fragmentation functions statistically. 



5. Conclusions 

Full jet reconstruction expands the kinematic reach 
to much larger values. These large momenta reaching 
50 GeV can be studied in heavy ion collisions for the 
first time. New physics effects should be considered 
when interpreting the results at large momentum. A 
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Figure 5: z dependence of the ratio of the fragmentation 
functions for recoil jets (normalized to the number of re- 
coil jets) measured in 0-20% central High Tower triggered 
Au+Au events to High Tower triggered p+p collisions [l5j]. 
The systematic uncertainties in the estimation of the width 
of the Gaussian parameterization of the background fluc- 
tuations are presented as solid envelopes. 



possible effect might be the momentum dependence 
of the relative contributions of quark and gluon sub- 
processes to inclusive jet production. In elementary 
p+p collisions these contributions vary with respect 
to the jet momentum (32j. The relative contributions 
might be even different in a heavy ion environment 
when a quark gluon plasma is produced, affecting the 
expected shape of the jet spectra and therefore of the 
nuclear modification factors. Another possibility is 
that at large momentum fraction x, initial state ef- 
fects (such as the EMC effect which is the deviation 
between structure functions of heavy ions to light ions) 
are observed to be as large as 15% [33]. There might 
be other effects like the EMC effect playing a major 
role in the relative suppression or enhancement of nu- 
clear modification factors at large momentum. 

The new Monte-Carlo based simulations of jet 
quenching in medium such as Jewel [34J, Q-Pythia 
|35j | and YaJEM [3™ and complementary analytic cal- 
culations [13, l38l. I39L [40| recently became available to 
pursue a quantitative analysis of jet quenching as ob- 
served in heavy ion collisions. However there are many 
uncertainties (e.g., how hadronization is treated) in 
the predictions of these models and calculations. To 
confront the calculations with data, new robust QCD 
jet observables that are unaffected by the pt cuts and 
hadronization need to be explored experimentally. For 
example the subjet observable is infrared safe and in- 
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sensitive to hadronization and may be used to study 
the jet quenching (34|. 

Unbiased reconstruction of jets in central heavy ion 
collisions at RHIC energies would be a breakthrough 
to investigate the properties of the hot QCD matter. 
The studies shown here indicate that reconstruction 
of jets with a uniquely large kinematic limit may in- 
deed be possible in heavy ion events. However jet 
reconstruction in heavy ion collisions is not yet free of 
biases. Biases are introduced due to selection of parti- 
cles such as pt cuts to reduce the fluctuations of heavy 
ion background, requirement of algorithmic parame- 
ters such as cone size or resolution parameter, and the 
collection of events with thresholds to enhance the jet 
rates. To assess fully the systematic uncertainties of 
jet measurements these effects must be investigated 
further. These results motivate us to look further into 
multiple channels for consistency checks (inclusive, di- 
jets, h-jets, gamma-jets to measure qualitatively new 
observables: ener gy f low, jet substructure, fragmenta- 
tion function [37l.l38l|.) 
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